Figure 1: Arithmetic means of crop speculative coverage.
coverage.lm1 <- lm(speculative.pct.coverage ~ as.factor(block) + species , data = coverage)
resid_panel(coverage.lm1)
joint_tests(coverage.lm1)
## model term df1 df2 F.ratio p.value
## block 3 33 8.267 0.0003
## species 11 33 13.630 <.0001
plot(emmeans(coverage.lm1 , "species"), comparisons = TRUE) +
xlab("Speculative coverage (%)" ) +
ylab("Species") #+ scale_y_discrete(labels = crops_labels)